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1.0  INTRODUCTION 


Detonation  physics  is  a  field  with  its  roots  in  the  investigations  of  combustion  wave 
properties  by  Hugoniot  (1887-1889),  Chapman  (1889)  and  others.1  A  solid  understanding  of  the 
development  and  consequences  of  these  and  other  fundamental  studies  is  crucial  to  any  research 
into  the  mechanics  of  detonation  and  into  the  subsequent  blast.  In  the  simplest  sense,  we  can 
regard  blast  as  an  interaction  of  the  detonating  explosive  with  its  surroundings,  e.g.,  the  air  or 
neighbouring  structures.  The  present  work  also  constitutes  a  brief  investigation  into  blast 
chemistry — the  interaction  of  detonation  products  with  the  air  and  with  themselves  to  form  other 
products.  The  detonation/blast  environment  is  highly  complex  due  to  the  presence  of  multiple 
species  and  chemical  processes  with  varying  rates.  A  proper  and  more  complete  understanding 
of  these  dynamics  may  lead  to  ways  of  controlling  blast  effects.  This  is  a  topic  of  significant 
interest  to  the  US  Air  Force.  The  present  work  begins  with  a  discussion  of  basic  detonation 
theory  to  orient  the  reader  to  our  topic  of  interest.  Moreover,  we  illustrate  that  conventional 
detonation  theory  treats  energy  release  as  a  single  lumped  chemical  parameter.  Although  we  will 
begin  to  depart  from  this  idea,  it  is  essential  to  understand  the  ideal  theory  in  that  it  provides  a 
simple  linkage  between  detonation  chemistry  and  blast  thermodynamics. 

A  computer  program  to  illustrate  this  theory  has  been  constructed  for  gaseous  detonation 
products  which  may  satisfy  either  of  two  equations  of  states.  The  first  corresponds  to  the 
calorically-perfect  gas  (CPG)  equation  of  state,  and  the  second  to  the  Jones-Wilkins-Lee  (JWL) 
equation  of  state.  The  first  is  chosen  because  it  is  a  simple,  yet  appropriate,  introduction  to 
solving  fundamental  detonation  physics  problems.  The  second  is  widely  used  and  has  been 
proven  reliable  when  practical  application  is  required,  so  incorporating  it  into  the  simple 
framework  previously  developed  results  in  a  comprehensive,  powerful,  and  simple  tool  for 
predicting  detonation  conditions.  This  hydrodynamic  analysis  makes  for  a  thorough  evaluation  of 
the  essential  physics  governing  the  detonation  problem  lending  immediately  to  practical 
application.  Model  validation  exercises  are  performed  using  five  well  known  explosive 
compounds.  Program  results  are  shown  to  correlate  well  with  experiments. 

As  a  second  technical  topic,  we  develop  a  subroutine  to  simulate  chemical  reactions  in 
the  post-detonation  environment.  This  subroutine  is  employed  by  LESLIE3D,  Large  Eddy 
Simulation  with  Linear  Eddy  modeling  in  3  Dimensions,  a  computer  program  developed  by 
Suresh  Menon  at  the  Georgia  Institute  of  Technology.'  It  simulates  the  time  evolution  of  the 
shock  wave,  keeping  up  with  such  flow  field  parameters  as  pressure,  temperature,  and  velocities 
while  performing  dynamic  large  eddy  simulation.  In  this  context,  large  eddy  simulation  (LES)  is 
a  family  of  techniques  that  divides  the  numerical  estimation  of  properties  into  two  parts.  First, 
filtered  governing  equations  are  solved  at  scales  larger  than  the  mesh  size.  The  solution  at  these 
“resolved”  scales  relies  upon  the  numerical  solution  of  the  governing  partial  differential 
equations.  Secondly,  scales  below  the  level  of  the  grid  cells  require  a  modeling  approach.  In  this 
case,  our  model  of  the  subgrid  stress  tensor  exploits  self-similarity  through  a  mathematical 
analogy  with  the  Leonard  stress.  Most  simulations  with  this  level  of  complexity  deal  only  with 
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the  fluid  products  of  detonation,  but  common  explosives  also  produce  soot,  which  is  unique  in 
that  it  is  a  solid  product.  It  is  hoped  that  further  understanding  of  the  formation  of  solid 
byproducts  may  provide  greater  insight  into  the  mechanisms  of  explosive  energy  loss.  In  this 
study,  soot  particles  are  considered  to  be  made  entirely  of  carbon.  Soot  develops  within  the  blast 
as  a  result  of  four  primary  stages:  nucleation,  surface  growth,  oxidation,  and  agglomeration.3  We 
have  already  tested  a  nucleation  model  for  soot,  but  observation  has  shown  that  young  soot 
particles  display  very  rapid  mass  growth  in  the  presence  of  soot  precursors  like  acetylene.3  This 
fact  motivates  the  present  study.  Here,  the  soot  formation  mechanism  is  enhanced  by  adding  a 
second  step  to  the  reaction  mechanism — surface  growth  due  to  the  adsorption  of  acetylene  on  the 
surface  of  the  particles.  It  takes  into  account  the  growing  diameter  size,  its  effect  on  the  surface 
growth  reaction  rate,  and  the  effect  that  aging  throughout  the  flame  reduces  the  reactivity  of  the 
soot  particles. 

Results  of  this  investigation  are  compared,  in  the  case  of  detonation  conditions,  to  well 
established  detonation  criterion  for  multiple  explosive  compounds,  and  in  the  case  of  the  soot 
formation  and  kinetics  along  a  shock  wave,  to  an  earlier  study  which  excludes  soot  surface 
growth.  Primarily,  it  is  hoped  that  future  studies  will  allow  new  explosives  to  function  with 
increased  efficiency,  particularly  by  maximizing  energy  release  per  unit  explosive  mass.  With 
this  motivation  in  mind,  adding  the  third  step  in  our  soot  model,  oxidation,  to  the  simulation  may 
increase  our  understanding  of  the  chemistry  required  to  tailor  the  loading  of  oxidizer  versus 
explosive  compound  in  the  material.  In  time,  advanced  simulation  techniques  like  those 
described  in  part  below  may  provide  insight  for  improving  the  manufacturing  and  loading  of 
metalized  explosive  components  such  as  aluminum  particles. 
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2.0  DETONATION  PHYSICS 


2.1  Technical  Approach 

For  a  one-dimensional  steady  flow  combustion  wave,  treating  the  products  and  reactants 
as  compressible  fluids  and  fixing  the  coordinate  system  to  the  detonation  wave  front,  the  basic 
conservation  equations  of  mass,  momentum,  and  energy  are  given  by 


Pouo  —  P 1U1 
Po  +  Po^o  =  Pi  +  pxu\ 


fro  +  y  =  fri  +  y 


(1) 

(2) 

(3) 


where  i  —  0,1  will  hereafter  represent  the  initial  and  final  states,  respectively;  pL  is  the  density 
of  the  mixture;  iq  is  the  velocity  of  the  mixture  with  respect  to  the  combustion  wave;  p,  is  the 
pressure  of  the  mixture,  and  /q  is  the  enthalpy  per  unit  mass  of  the  mixture.  These  three 
equations  are  sufficient  to  define  the  Rayleigh  Line  , 


Pl-Po 
Vi-  v0 


2  2  2  2 
=  -p  o“o  =  “PlX 


(4) 


where  iq  =  1  /pt  is  the  specific  volume.  Equations  (1)  through  (4)  may  also  be  used  to  define 
the  Hugoniot  Equation  ,  a  relationship  completely  independent  of  the  equation  of  state,  i.e., 

fri  “  fro  -  q  =  \(Pi  ~  Po)Oo  +  ^i)  (5) 

where  q  —  h^Q  —  hf  is  the  difference  between  the  enthalpies  of  formation  of  reactants  and 
products.  When  q  is  positive,  (5)  is  referred  to  as  the  detonation  Hugoniot.  When  q  =  0,  this 
equation  relates  thermodynamic  states  in  unreacted  material  and  is  denoted  the  shock  Hugoniot. 
As  an  example,  the  shock  pressure  for  a  calorically  perfect  gas  may  be  written  as 


Pi  = 


(r  rro+i'j  _vi  +  zqi" 

IL'o-il  l'o  1  p  0 

[(“as-1!  r 


(6) 


where  Yi  is  given  by 


Yi  = 


(7) 


R  =  8.314  J/(mol  ■  K)  is  the  universal  gas  constant,  and  Cp  j  is  the  specific  heat  capacity  of  the 
mixture  in  J/ (mol  ■  K).  It  is  written  as 


Cp,i  =  ZiXyC 


u  p.y 


(8) 
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Figure  1:  Shock  Hugoniot  for  the  Combustion  of  Acetylene  using  the  Calorically  Perfect  Gas  Equation  of  State 


Shock  Hugoniot 


where  the  summation  runs  over  the  j  species  for  reactants  (i  =  0)  and  products  (i  =  1);  C p  is 
the  specific  heat  capacity,  Cp,  in  J/(mol  ■  K)  of  the  jth  species;  and  Xtj  —  Pij/pi  is  the  mass 
fraction  of  the  jth  species,  where  pLJ  is  the  partial  density  of  that  species.  The  simple  case  of  the 
combustion  of  acetylene,  shown  below,  is  considered. 

C2H2  +^02  2C02  +  H20  (9) 

For  the  combustion  of  acetylene,  a  graph  of  equation  (6)  produces  a  shock  Hugoniot  as  seen  in 
Figure  1.  With  qc2H2  —  1255.5  kj/mol4,  the  detonation  Hugoniot  for  this  reaction  is  produced 
and  shown  with  the  shock  Hugoniot  in  Figure  2. 

Now  that  the  basic  procedure  is  established,  we  may  develop  a  Hugoniot  for  a  more 
sophisticated  equation  of  state.  The  JWL  equation  of  state  may  be  written  as  5 

p  =  1 4  (l  -  ^)  exp( -RtV)  +  B  (l  -  exp (~R2V)  +  y  (10) 


where  A,B,w,  and  R1/2  are  unique  constants  specified  for  each  explosive,  V  —  v1/vQ  ,  and  e  is 
the  volumetric  internal  energy  and  is  given  by  e  —  p0E,  where  E  is  specific  internal  energy. 
Now  the  form  of  equation  (5)  must  include  specific  internal  energy  instead  of  sensible  enthalpy 


Ei  ~  E0  -  q  =  \  (pi  +  p0)(v0  -  vx) 


Solving  (10)  for  E  and  substituting  into  (11),  we  obtain6 


PlOO  = 


{ l  [  4(l-fiw)exp(-fiil/)+B(1~^T)exp(-fi2t')  M  ^ ^]+ 1 


(ii) 


(12) 
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Figure  2:  Detonation  and  Shock  Hugoniots  for  the  Combustion  of  Acetylene  using  the  Calorically  Perfect  Gas  Equation  of 
State 

where  E0  represents  the  internal  energy  of  the  solid  explosive,  and  is  given  by 

E0  =  CVT0  (13) 

where  T0  is  the  initial  temperature,  and  Cv  is  its  constant  volume  specific  heat  capacity. 
Approximate  values  of  Cv  are  used  for  the  explosives  of  interest.  Given  an  initial  temperature  of 
about  103K  and  a  specific  heat  capacity  for  an  average  explosive  of  about  103J/(kg  ■  K)  ,  E0  is 
on  the  order  of  105J/kg.  Compared  to  the  other  terms  in  (12),  its  contribution  to  pressure  is 
negligible,  so  E0  may  be  set  to  zero.  Figures  3  through  7  show  the  JWL  detonation  Hugoniots 
for  Trinitrotoluene  (TNT),  Cyclo-tetramethylene  tetranitramine  (HMX),  a  Cyclo-trimethylene 
trinitramine  or  RDX  derivative  (Composition  C-4),  a  Triaminotrinitrobenzene  derivative  (PBX- 
9502)  and  Pentaerythritol  tetranirate  (PETN),  respectively.  The  JWL  equation  of  state 
parameters  are  provided  in  Table  l.7 


Table  1.  JWL  Equation  of  State  Parameters 


Explosive 

A  (Mbar) 

B  (Mbar) 

Ri 

r2 

CD 

TNT 

3.712 

0.3231 

4.15 

0.95 

0.30 

HMX 

7.783 

0.07071 

4.20 

1.00 

0.30 

C-4 

6.0997 

0.1295 

4.5 

1.4 

0.25 

PBX-9502 

4.603 

0.9544 

4.0 

1.70 

0.48 

PETN 

5.731 

0.20160 

6.00 

1.80 

0.28 

In  the  interest  of  comparison,  a  products -based  shock  Hugoniot  is  shown  for  each 
explosive;  its  locus  is  obtained  by  setting  q  —  0  in  the  product  JWL  Hugoniot  formula  (12).  In 
doing  so,  we  can  illustrate  the  shift  in  the  Hugoniot  locus  caused  by  the  detonation  energy  term. 
In  practice,  solid  explosive  materials  require  different  equations  of  state  (usually  not  the  JWL 
equation)  in  order  to  capture  the  correct  response  to  shock  pressure.  This  solid  phase  equation 
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Detonation  and  Shock  Hugoniots,  TNT 


Figure  3:  Detonation  and  Shock  Hugoniots  for  TNT  using  the  JWL  Equation  of  State 
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Figure  4:  Detonation  and  Shock  Hugoniots  for  HMX  using  the  JWL  Equation  of  State 


Detonation  and  Shock  Hugoniots#  W1X 
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Detonation  and  Shock  Hugoniots,  Conp  04 


Figure  5:  Detonation  and  Shock  Hugoniots  for  Composition  C-4  using  the  JWL  Equation  of  State 


Detonation  and  Shock  Hugoniots,  PBX-9582 


Figure  6:  Detonation  and  Shock  Hugoniots  for  PBX-9502  using  the  JWL  Equation  of  State 
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Detonation  and  Shock  Hugoniots,  PETN 


Figure  7:  Detonation  and  Shock  Hugoniots  for  PETN  using  the  JWL  Equation  of  State 

(e.g.,  a  Mie-Griineisen  or  Hayes  formula,  suitably  calibrated  for  the  explosive)  is  substituted  in 
(12)  with  q  =  0  to  obtain  the  shock  Hugoniot  for  the  solid  explosive. 

With  the  background  provided  on  the  conservation  laws  and  on  the  equation  of  state  for 
the  detonation  products,  we  may  now  determine  the  properties  of  the  flow  field  at  detonation 
conditions.  In  the  ideal  theory,  detonation  is  modeled  as  an  instantaneous  phase  transition.  The 
solid  explosive  originally  at  the  state  P0,  pQ  transitions  directly  to  the  steady  state  detonation 

conditions  with  pressure  and  density,  Pi  and  p] ,  respectively.  That  is  to  say,  the  phase 
transformation  does  not  move  along  the  Rayleigh  line  in  time.  The  change  in  properties  is 
effectively  discontinuous.  Any  realizable  detonation  condition  must  lie  on  both  the  Detonation 
Hugoniot  and  the  Rayleigh  Line  as  shown  in  Figure  8.  The  two  solutions  shown  correspond  to 
strong  and  weak  detonations;  however,  the  strong  detonation  state  is  rarely  observed.  The  weak 
detonation  is  generally  not  observed.1  The  most  commonly  observed  detonation  condition  is 
Chapman-Jouguet  (CJ)  detonation.  This  state  exists  where  the  Rayleigh  Line  is  tangent  to  the 
Detonation  Hugoniot,  a  position  referred  to  as  the  CJ  Point.4  It  is  also  the  state  found  in  steady 
detonations.  The  existence  of  the  CJ  point  is  not  deduced  from  the  conservation  laws.  Rather  its 
existence  is  explained  from  a  conjecture  involving  detonation  physics  substantiated  by  empirical 
measurements.1  This  point  may  be  found  by  plotting  a  function  F,  defined  as  the  difference 
between  the  Rayleigh  line  and  detonation  Hugoniot  slopes,  and  finding  its  zero.  A  computer 
program  is  used  to  split  up  the  plot  into  regions  AV  wide  in  order  to  determine  at  which  V  the 
sign  of  F  changes,  indicating  a  root.  The  function  F  is  determined  to  be 

r  _  PiQO-Po  Pi{v+)-Pi(-V~ )  (U) 

V-l  AV  '  ' 
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Detonation  Huconiot  and  Ra«jlei|(h  Line 


Figure  8:  Two  Possible  Detonation  States  Identified  by  Intersection  of  Rayleigh  Line  and  Detonation 
Hugoniot 

where  Pt(V )  is  given  by  (12).  In  equation  (14),  we  have  introduced  some  new  notation.  Let 

Pi<y±)  =  pi(y±f),  (is) 

where  AV  is  given  by 

AV  =  (16) 

k 

Vmax  is  a  reasonable  upper  bound  on  the  value  of  V ;  Vmin  is  a  reasonable  lower  bound  on  the 
value  of  V,  and  k  is  the  number  of  sections  comprising  the  region  between  Vmax  and  Lmin.  The 
parameter  k  is  assigned  a  value  of  100,000  in  order  to  accurately  the  value  of  VCJ  where 

F(VCJ )  =  0  as  in  (14).  In  the  next  section,  this  algorithm  is  employed  to  calculate  the  Chapman- 
Jouguet  detonation  conditions  for  a  series  of  explosives. 


2.2  Results 

The  simple  root  finding  algorithm  (14-16)  developed  in  the  preceding  section  has  been 
used  to  determine  the  CJ  conditions  for  the  legacy  explosives  TNT,  HMX,  Composition  C-4, 
PBX-9502  and  PETN.  The  CJ  points  are  graphed  on  the  detonation  Hugoniots  for  these  materials 
in  Figures  9  through  13,  respectively,  for  the  subject  explosives.  The  CJ  pressures  are  validated 
exhibiting  excellent  agreement  with  archival  data.  This  information  is  provided  in  Table  2. 
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CJ  Condition,  TNT 


Figure  9:  Chapman-Jouguet  State  for  TNT 


Figure  10:  Chapman-Jouguet  State  for  HMX 
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GJ  Condition,  Conn  C-4 


Figure  11:  Chapman-Jouguet  State  for  Composition  C-4 


Detonation  and  Shock  Hufponiots,  P8X-9502 


Figure  12:  Chapman-Jouguet  State  for  PBX-9502 
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CJ  Condition,  PETM 


Figure  13:  Chapman-Jouguet  State  for  PETN 


Table  2.  Validation  of  Numerically  Calculated  Chapman-Jouguet  Pressures 


Explosive 

CJ  Pressure  (GPa) 

Calculated 

Reported7 

TNT 

21.0 

21.0 

HMX 

42.2 

42.0 

Composition  C-4 

28.0 

28.0 

PBX-9502 

30.2 

30.2 

PETN 

14.0 

14.0 

It  is  evident  that  this  algorithm  generates  CJ  pressures  matching  the  known  values  for  the 
subject  explosives  within  one-half  percent  error.  The  largest  discrepancy  is  associated  with  the 
CJ  pressure  for  HMX.  Our  prediction  is  42.0  GPa,  while  the  value  reported  by  the  previous  study 
is  42.2  GPa.  Dobratz  reports  that  the  archival  CJ  pressure  for  HMX  is  estimated  ,  so  the 
disagreement  between  this  value  and  the  calculated  pressure  is  minor.  This  study  is  important 
because  of  its  fundamental  nature  and  its  wide  range  of  applicability.  The  Hugoniot  equation, 
Rayleigh  Line,  and  CJ  conditions  all  arise  naturally  from  the  conservation  laws  and  the  JWL 
equation  of  state.  Although  the  JWL  equation  of  state  remains  controversial,  it  still  remains  as  a 
workhorse  in  performing  this  type  of  calculation,  and  it  is  a  natural  extension  of  the  ideal  theory 
based  upon  a  single  “measurement”  of  detonation  energy  (unlike  the  combustion  of  acetylene). 
In  the  next  section,  we  introduce  finite  rate  blast  chemistry  where  each  chemical  reaction 
existing  in  the  blast  has  its  own  specification  for  reaction  energy. 
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3.0  SIMULATING  DETONATION  SOOT  FORMATION 


3.1  Technical  Approach 

The  preceding  section  presents  ideal  detonation  theory  from  the  standpoint  of  combustion 
because  detonation,  like  combustion,  is  a  complex  chemistry  problem.  There  are  two  major 
differences  between  combustion  and  detonation.  In  the  first  place,  the  chemical  rates  for 
detonation  are  far  more  rapid  than  those  for  combustion.  In  fact,  at  the  laboratory  scale, 
detonation  reactions  appear  to  be  nearly  instantaneous.  Secondly,  combustion  processes  tend  to 
acquire  oxygen  from  their  surrounding  environment.  Explosives  contain  an  amount  of  oxidizer 
sufficient  to  conduct  the  detonation.  When  “burned”,  both  fuels  and  explosives  produce  a  wide 
array  of  products.  Even  simple  fuels  like  methane  have  reaction  mechanisms  involving  eight  or 
more  species.  More  complex  molecules  such  as  HMX  may  involve  one  hundred  products  or 
more.  These  products  include  both  gases  and  solids.  The  family  of  solid  products  denoted  as 
“soot”  has  acquired  a  certain  “enigmatic”  nature. 

Neither  the  composition  of  nor  the  mechanism  behind  the  formation  of  soot  via 
combustion  is  well  understood.  For  our  purposes,  soot  is  assumed  to  be  purely  solid  carbon,  but 
its  formation  seems  to  rely  on  the  presence  of  one  or  more  chemical  precursors,  e.g.,  acetylene. 
This  study  is  based  on  a  four  step  model:  nucleation,  surface  growth,  oxidation,  and 
agglomeration.  The  first  two  steps  govern  the  mass  formation,  and  the  last  two  deal  with  the 
interaction  of  the  soot  after  formation.  The  present  investigation  concerns  itself  with  the  mass 
formation  of  the  soot — the  first  two  steps.  All  calculations  are  done  for  each  volume  element  of 
the  simulation  grid  (or  mesh),  so  many  of  the  units  given  below — densities,  concentrations, 
etc. — are  for  mixture-based  quantities. 

Some  studies  suggest  that  of  particular  importance  to  soot  growth  is  the  presence  of 
acetylene,  C2H2,  which  correlates  with  a  very  rapid  initial  growth  phased  In  this  case,  the  soot 
precursor,  acetylene,  is  assumed  to  arise  from  Lauryl  methacrylate  (LMA),  a  component  of  the 
HMX-based  explosive.  A  further  assumption  made  is  that  some  active  soot  (carbon)  nuclei  are 
formed  from  the  breakdown  of  acetylene 

C2H2  ->  2C(s)  +  H2  (17) 

The  reaction  rate,  then,  is  given  as3 

I?i  = /Ci(T)[C2H2]  (18) 

where  is  given  in  mol/ (cm3  ■  s),  and  [C2H2]  is  the  current  concentration  of  acetylene  in  the 
mixture  and  is  given  in  mol/cm3.  The  rate  constant  k1  (71),  in  s-1,  is  given  by3 

kt(T)  =  Aexp(Ea/RT)  (19) 

where  Ea  is  the  activation  energy  and  R  is  the  gas  constant  for  acetylene.  The  literature  states  has 
Ea/R  =  21,100  K  and  the  pre-exponential  factor  A  —  10,000  s-1  for  a  C2H4  (ethylene)  flame 
burning  with  an  oxidant  stream  of  22%  02  and  78%  N2.  Note  that  ethylene  is  a  chemical 
precursor  for  acetylene.  The  second  step  is  the  adsorption  of  C2H2  on  the  surface  of  the  soot 
particles3,  i.e., 
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C2H2  +  nC(s)  ->  (n  +  2)C(s)  +  H2 


(20) 


An  arbitrary  soot  particle  size  n  is  chosen  to  begin  this  step  of  the  soot  formation  process.  An 

o 

initial  particle  size  has  been  reported  as  20  to  50  carbon  atoms,  so  n  is  chosen  as  20.  This 
reaction  is  again  governed  by  the  acetylene  concentration,  but  also  by  the  number  of  active 
adsorption  sites,  which  is  related  to  the  surface  area  of  the  particle.  Studies  have  shown  that,  on 
the  average,  older  particles  are  less  reactive  than  younger  particles. 1  To  account  for  this  decrease 
in  reactivity  with  age,  the  number  of  reaction  sites  is  assumed  to  be  proportional  to  the  square 
root  of  the  surface  area  available  locally  in  the  flame.  So  the  governing  reaction  rate  is 


R2  =  fe2(n/(5)[C2H2]  (21) 

o 

where  [C2H2]  is  the  current  concentration  of  acetylene  in  the  mixture  and  is  given  in  mol/cm  ; 
/(S)  is  the  square  root  of  the  available  surface  area  S,  and  k2(T )  is  the  rate  constant.  Based 

2  3  3 

again  on  the  mixture,  S  is  given  in  units  of  cnr/cm  and  is  written  as 


S  =  7t(dp)(plV) 


(22) 


where  p,  in  g/cnr  ,  is  the  mixture  density.  Also,  N  is  the  total  number  of  soot  particles  per  gram 
of  mixture,  and  dp  is  the  diameter  of  each  soot  particle,  given  as 


_  (  6  1  YC(s)  \ 
V  71  Pc  N  ) 


(23) 


Note  that  pc  =  2  g/ cm3  is  the  soot  density  (based  upon  a  reasonable  value  for  the  density  of 
carbon  ),  and  YC(S)  is  the  mass  fraction  of  solid  carbon  in  the  mixture.  To  create  a  more  realistic 
soot  model,  it  is  desirable  to  build  a  number  density  evolution  equation  into  LESLIE3D  as  a  part 
of  its  governing  equations.  This  equation  actually  captures  the  evolution  of  the  number  of  soot 
particles  per  unit  volume  in  time  and  space.  Implementing  this  method  requires  further 
development  of  LESLIE3D  and  is  reserved  for  future  research.  At  present,  we  assume  that  N  is 
fixed  at  100  particles/g  of  mixture.  The  mass  fraction  is  given  by 


Yc(s)  = 


[C(s)]Mc(s} 

P 


(24) 


where  [C(s)]  is  the  concentration  of  solid  carbon  in  mol/cm3,  and  MC(S)  is  the  molar  mass  of 
solid  carbon,  12.011  g/mol  .  From  equations  (22-24),  equation  (21)  becomes^ 


R2  =  /c2(T)[C2H2](7i)?(^[C(s)])i(piV)5  (25) 

where  the  rate  constant  k2(T ),  in  cml/2/cm/s  for  the  soot ,  is  given  by3 

k2(T)  =  Aexp(-Ea  /  RT)  (26) 

where  Ea/R  =  12,100  K  and  the  preexponential  factor  A  =  6,000  cm  /cm/s  has  been  determined 
for  a  C2H4  flame  burning  with  an  oxidant  stream  of  22%  02  and  78%  N2.  Two  other  reactions, 
the  combustion  of  carbon  monoxide9  and  the  combustion  of  hydrogen10,  which  do  not  directly 
produce  soot  are  included  in  this  model,  i.e., 
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2CO  +  02  ->  2C02 

(27) 

2H2  -I-  02  ->  2H20 

(28) 

For  these  reactions,  the  pre-exponential  factors  (A)  are  3.98  x  1014  s_1 
Ea/R  is  set  equal  to  20,130  K  and  17,614  K,  respectively.9'10  The  entire 
this  model  may  be  written  as 

and  1.8  x  1013  s_1. 
mechanism  used  for 

C2H2  ^2C(s)  +  H2 

(28) 

C2H2  +  nC(s)  ->  (n  +  2)C(s)  +  H2 

(29) 

2CO  +  02  ^  2C02 

(30) 

2H2  +  02  ->  2H20 

(31) 

with  reaction  rate  information  given  in  Table  3. 3  910 

Table  3.  Reaction  Rate  Information 


Reaction 

A 

Ea/R  (K) 

Soot  Nucleation 

10,000  s'1 

21,100 

Soot  Surface  Growth 

6,000  cm3/2/cm/s 

12,100 

CO  Combustion 

3.98xl014  s1 

20,130 

H2  Combustion 

1.8xl013  s'1 

17,614 

A  subroutine  representing  this  reaction  mechanism  has  been  written  in  standard  CHEMKIN 

CKWYP  form  and  integrated  into  LESLIE3D.11  LESLIE3D  calls  the  CKWYP  subroutine  and 

12 

solves  for  the  species  concentrations  using  its  locally  dynamic  subgrid  kinetic  energy  model. 
The  filtered  governing  equations  (including  the  species  equations)  are  closed  with  the  use  of  the 
Eddy  Break-Up  turbulent  chemical  reaction  closure  model.  This  model  ensures  that  the 
reaction  rates  are  properly  limited  for  turbulent  mixing.  Without  this  type  of  closure,  the  kinetic 
reaction  rates  are  over-predicted.  The  following  section  of  the  report  contains  the  results  of  this 
model  applied  to  a  typical  explosive. 


3.2  Results 

The  reaction  mechanism  described  in  Section  3.1  has  been  applied  to  the  blast  field 
produced  by  the  detonation  of  an  HMX-based  explosive.  Although  LESLIE3D  cannot 
accomplish  the  actual  detonation,  the  detonation  is  conducted  in  a  separate  computer  code  and 
imported  into  LESLIE3D.  LESLIE3D’s  simulation  begins  at  8.5  ps  of  problem  time.  We 
compare  the  results  obtained  for  the  pure  nucleation  model  with  those  of  the  latest  model  that 
includes  nucleation  with  the  surface  growth  reaction  step.  Figure  14  shows  the  total  number  of 
carbon  atoms  in  the  simulation  with  respect  to  time  for  pure  nucleation,  while  Figure  15  shows 
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Carbon  Fornation  with  Tine,  Nucleation 


Tine  (seconds) 


Figure  14:  Soot  Production  via  Pure  Nucleation  for  an  HMX-based  Explosive 


Carbon  Fornation  with  Tine,  Surface  Growth 


Tine  (seconds) 


Figure  15:  Soot  Production  via  Nucleation  Plus  Surface  Growth  for  an  HMX-based  Explosive 

the  total  number  of  carbon  atoms  in  the  simulation  with  respect  to  time  for  the  inclusion  of  terms 
representing  surface  growth  through  adsorption.  As  we  may  expect  from  an  examination  of  the 
reaction  mechanism,  the  soot  nucleation  process  is  quite  strong  in  the  early  moments  of  the  blast. 
(Note  that  this  model  contains  no  soot  destruction  mechanism).  Yet  the  slopes  of  the  curves  in 
Figures  14  and  15  begin  to  reduce  at  about  15  ps.  This  effect  is  likely  due  to  the  relatively  small 
amount  of  acetylene  used  to  initialize  the  simulation.  Recall  that  our  original  acetylene 
concentration  is  based  upon  an  assumed  set  of  decomposition  products  for  Lauryl  methacrylate,  a 
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Carbon  Fornation  with  Tine 


Tine  (seconds) 

Figure  16:  A  Comparison  of  Soot  Production  via  Nucleation  and  Nucleation  Plus  Surface  Growth  for  an 
HMX-based  Explosive 


Carbon  Fornation  with  Tine 


Tine  (seconds) 

Figure  17:  A  Comparison  of  Soot  Production  via  Nucleation  and  Nucleation  Plus  Surface  Growth  for  an 
HMX-based  Explosive  over  the  Time  Period  of  16  ps  to  17  ps 

constituent  of  the  explosive.  Figure  16  shows  both  plots  for  pure  nucleation  and  surface  growth 
in  order  to  facilitate  direct  comparison.  The  difference  between  these  two  curves  is  small,  but  the 
surface  growth  term  does  increase  soot  production,  especially  at  later  times.  By  magnifying  the 
plot,  may  obtain  a  clearer  view  of  the  differences.  Figure  17  shows  a  magnified  view  of  soot 
production  from  16  to  17  ps  while  Figure  17  shows  soot  production  over  the  time  period  of  25 
to  27.5  ps.  At  16.5  ps,  the  difference  (for  the  entire  blast  field)  is  approximately  2.4  g  of  carbon. 
However,  the  increase  in  surface  growth  causes  the  production  of  carbon  to  escalate  to  about  14 
g  at  26.5  ps.  Bearing  the  number  of  assumptions  in  mind,  it  is  interesting  to  see  the  quantitative 
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Carbon  Fornation  with  Tine 


Tine  (seconds) 


Figure  18:  A  Comparison  of  Soot  Production  via  Nucleation  and  Nucleation  Plus  Surface  Growth  for  an 
HMX-based  Explosive  over  the  Time  Period  of  25  ps  to  27.5  ps 

differences  in  soot  production  predicted  by  this  simple  model.  In  a  follow-on  step,  it  would  be 
instructive  to  see  the  effects  of  including  soot  oxidation  within  the  reaction  mechanism.  This  process  is 
relatively  easy  to  integrate  into  the  overall  reaction  mechanism.  Yet,  it  is  perhaps  of  greater  importance  to 
first  improve  the  selection  of  detonation  products.  The  present  set  of  detonation  products  is  too  simple. 
Many  possible  byproducts  of  HMX’s  detonation  are  overlooked.  The  presence  of  more  products  may 
significantly  alter  blast  chemistry. 
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4.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  present  work  takes  some  preliminary  steps  in  considering  the  processes  of  detonation  and  the 
subsequent  burn  of  detonation  products  in  the  blast  field  in  a  unified  way.  Detonation  has  most  often  been 
treated  from  the  standpoint  of  shock  physics.  The  detonation  wave  is  often  considered  to  be  a 
discontinuity  propagating  through  the  explosive.  We  have  implemented  this  theory  and  have  shown  that  it 
produces  good  results  for  pressure.  Yet,  it  is  important  to  realize  that  detonation  is  a  chemical  reaction 
that  takes  place  over  a  very  small  region  of  space.  This  reaction  creates  products,  chemical  species  that 
react  with  the  surroundings.  Moreover,  these  reactions  create  solid  carbonaceous  products  collectively 
called  “soot”.  The  second  part  of  this  study  considers  a  proposed  mechanism  for  the  production  of  soot  in 
the  blast  field  for  an  HMX -based  explosive.  This  mechanism  is  tested  using  the  LESLIE3D  multiphase 
physics  computer  code.  Using  turbulent  chemistry,  LESLIE3D  produces  time  accurate  curves  for  soot 
production  starting  from  an  initial  detonation  solution.  The  model  gives  good  qualitative  results — more 
soot  is  being  produced  when  the  surface  growth  step  is  included,  and  we  observe  a  proper  limiting  of 
carbon  production  as  the  initial  carbon  field  begins  to  disperse  and  rarefy.  The  difference  in  the  actual 
number  of  particles  between  the  two  models  seems  relatively  small,  but  does  result  in  an  appreciable 
increase  in  soot  mass.  Remember  that  even  in  regular  combustion  processes,  soot  production  is  often  low. 
It  is  often  detected  as  powdery  smears  on  engine  exhaust  piping.  It  is  also  apparent  that  fixing  the  total 
number  of  soot  particles  at  100  makes  for  a  low  estimate,  as  the  actual  number  of  particles  is  on  the  order 
of  1025  for  most  of  the  simulation.  Moreover,  we  know  that,  in  reality,  this  number  actually  changes. 
Representing  this  aspect  of  the  physics  requires  a  better  soot  model  incorporating  an  evolution  equation 
for  the  soot  particle  number  density.  One  could  also  expand  the  overall  reaction  mechanism  for  detonation 
products  to  include  soot  oxidation,  soot  agglomeration,  and  governing  reactions  for  other  by-products  of 
the  HMX -based  explosive.  One  could  represent  soot  transport  by  including  massive  Lagrangian  particles. 
These  particles  can  propagate  through  the  flow  field  under  inertial  and  drag  forces.  This  aspect  of  the 
model  may  provide  better  ideas  of  soot  dispersion.  That  is  to  say,  soot  that  disperses  more  quickly  cools 
more  rapidly.  It  is  also  prudent  to  study  the  overall  extraction  of  energy  release  data  throughout  the  blast 
field  to  obtain  greater  insight  into  the  effect  of  these  reactions  on  the  overall  blast  chemistry.  The  best 
recommendation  for  the  enhancement  of  this  research  is  that  these  numerical  simulations  continue  in 
parallel  with  an  experimental  investigation  into  both  the  gaseous  detonation  products  as  well  as  into 
identification  of  the  soot’s  chemical  make-up  created  by  explosives  of  choice. 
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